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Recently, Papadopoulos et al., 2010 using sub-mm CO molecular line observations of nearby ultra- luminous IRAS 
galaxies, (U)LIRGs with Lir > 10 11 Lq, have found that exceptionally large gas column densities (Nh > 10 25 cm~ 2 ) 
can be present across some of the very dense gaseous disks that are typically found in these objects. They also proposed 
a diagnostic for finding such sources using CO and HCN molecular lines. Given that such high column densities are 
expected to absorb any X-ray luminous AGN, yielding Compton-thick sources, we set out to explore whether this can 
be discerned using X-ray observations. More specifically we examine X-ray spectral observations of 14 sources in their 
sample, using public Chandra observations (0.5-10 keV) for eleven sources as well as BeppoSAX results (2-100 keV) 
from the literature for another three sources. Our goal is to find candidate Compton-thick AGN and to check whether 
the molecular line selection criterion is successful in selecting such systems. X-ray spectroscopy reveals four candidate 
Compton-thick AGN of which half fall within the high obscuration region in the molecular line ratio diagnostics. Of the 



remaining five sources falling into the 'high dust obscuration' box, one (Mrk273) is highly obscured (Nh 



2 ) 



while in the other four the X-ray emission is most probably associated with star-forming processes rather than an AGN 
on the basis of their X-ray and mid-infrared properties. Overall, we argue that although this method as expected 
cannot recover all Compton-thick AGN, there are no examples of X-ray luminous AGN inside that region that have low 
obscuration, suggesting that this method is efficient in finding heavily obscured AGN in dust-enshrouded star-forming 
galaxies. The above results bear important implications for future joint ALMA and X-ray observations for the detection 
of Compton-thick AGN. 

Key words. X-rays: general; X-rays: diffuse background; X-rays: galaxies; Galaxies: active; Submillimeter: galaxies; 
Infrared: galaxies 
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1. Introduction 

X-ray sur veys provide the most effic ient tool for identify- 
ing AGN (Brandt & Hasinger, 2005). This is because the 
hard X-rays can penetrate large amounts of dust and gas 
and thus can reveal hidden AGN which would be very dif- 
ficult to detect at optical wavelengths. However, even the 
hardest X-rays cannot easily escape fro m very large absorb 
ing c olumn densities, > 10 24 cm" 



( Murphy fc Yaqoob 
2009 1 . A column density of a few times 10 z4 cm _z absorbs 
X-rays up to 20 keV, while for column densities higher than 
. — 10 25 cm -2 practically all X-ray photons are attenuated. 
Such immense column densities (> 10 24 cm~ 2 ) correspond 
5_i to the so-called Compton-thick AGN, where Compton seat- 
ed tering becomes the dominant attenuation process. In these 
cases, reflected emission from the backside of the obscuring 
screen, in a form of a torus, remains the only evidence of 



the AGN presence in X-rays (Matt et al. 2004 



:4atn 



Such extreme column densities could be relatively eas- 
ily explored in the IR or sub-mm part of the spectrum, 
through their thermal or molecular emission. For exam- 
ple, sub-mm observations have drawn much attention re- 
centl y as a means of discov ering heavily obscured AGN 
(e.g. Alexander et al. 2005). The discovery of significant 



dust optical depth s even at submillimete r wavelengths (in 
the case of Arp220; Sakamoto et al. 2008 1 reveals that such 
highly obscured systems ma y not be uncommon. Recently, 
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Papadopoulos et al. (2010) (hereafter P10) made the in- 
teresting suggestion that observations of suitable molec- 
ular lines can reveal cases of very high column densities 
(and thus possible obscuration) present in some reservoirs 
of molecular gas in LIRGs. The average molecular gas den- 
sities in some of these systems can be so high that if the 
gas is distributed in a gaseous disk it can yield extremely 
high column densities and dust optical depths that sup- 
press the high-J (6-5) CO line emission to levels well be- 
low those expected in starburst galaxies such as ULIRGs. 
In such cases, if an active galactic nucleus were embedded 
in the (U)LIRG, the line of sight towards it would also 
be severely obscured making this source Compton-thick. 
A molecular line diagnostic to identify such severely ob- 
scured sources was proposed by P10. This helps distinguish 
between the degenerate possibilities of high dust obscura- 
tion and large reservoirs of non star-forming cool and dif- 
fuse gas, whose high-J (6-5) CO lines are genuinely faint 
rather than immersed to a strong dust continuum. They 
suggest that high fractions of very dense molecular gas 
(with n(H2) > 10 4 cm~ 3 ) can be traced by a high HCN(1- 
0)/CO(1-0) ratio (R>0.1). Then a low (R < 0.3) CO(6- 
5)/(3-2) ratio can only be the result of high dust optical 
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depths (r434 Mm <~ 1). Consequently there should be a very 
high probability of Compton-thick AGN for (U)LIRG hosts 
lying inside their aforementioned selection "window" . On 
the other hand, we note that gaseous disks (in the center of 
which the AGN are presumably embedded) with low aver- 
age column densities can often yield large column densities 
along the particular line of sight towards the AGN as a re 



suit of the i r high ly turbulent nature (Wada, Papadopoulos 
& Spaans 20091. Thus Compton-thick sources can be ex- 



pected also outside the dust-obscuration selection window, 
the vast majority of sources inside it however should be 
Compton thick. In other words, if molecular gas disks have 
on average high column densities, as the supression of CO 
J=6-5 emission from their star-forming gas would imply, 
this will almost certainly be the case for the special line of 
sight towards the AGN as well, but the reverse does not 
necessarily hold. Thus the selection window of highly dust 
obscured star-forming systems advocated by P10 ought to 
select almost exclusively Compton-thick AGN, but not all 
of them. 

In this paper we are examining the X-ray properties of 
the systems in P10, using both literature spectra as well as 
Chandra archival observations. Our goal is to test whether 
AGN residing in the P10 box have a higher probability 
for being classified as Compton-thick sources in X-rays. 
We adopt H D = 75kms- 1 Mpc" 1 , M = 0.3, fi A = 0.7 
throughout the paper. 

2. Data & analysis 

2.1. The sample 

P10 present molecular line spectra for 20 (U)LIRGs. 
For six objects there are no HCN observations avail- 
able so that the HCN/CO vs. CO(6-5)/CO(3-2) diag- 
nostic cannot be applied. For the remaining 14 objects 
there are X-ray spectra available either in the literature 
or in publicly available Chandra observations. Although 
in many cases the Chandra observations may have been 
presented elsewhere, we chose to re-analyze the data 
in order to provide a uniform analysis of the sample. 
To our knowledge, in six cases X-ray spectral fits to 
the Chandra data are presented here for the first time: 
VII Zw 031, IRAS 10565+2448, IRAS 12112+0305, Arp 193, 
IRAS 17208-0014, IRAS 23365+3604. However, X-ray lumi- 
nosities or hardness ratios derived from Chandra observa- 
tions have been presented before in the literature for some 
of those (see Tablejl]). In the cases of Mrk231, NGC6240 
and NGC 7469, where there are high signal-to-noise spec- 
tra in the 2 — 100 keV band reported in the literature from 
BeppoSAX, we prefer to present that data instead, as high 
energy observations can be much more efficient in detecting 
Compton-thick sources. 



2.2. Spectral fits 

We fit the Ch andra data in the 0.5-10 keV band using 
XSPECv. 12.6 ( |Arnaud|[T996| . Initially, we fit an absorbed 
power-law model to the spectra. For more complex spec- 
tra, and judging on the basis of the data residuals from 
the model fit, we additionally fit a thermal component 
(Raymond-Smith) to account for X-ray emission from star- 
forming processes, and/or an FeKa line around 6.4 keV. 
The C-statistic has been employed (Cash 19791. The er- 



rors correspond to the 90% confidence level as is custom- 
ary in X-ray Astronomy. The X-ray spectral properties and 
luminosities are presented in Table[T] In the Chandra data 
presented here, an FeKa line has been detected in the case 
of the sources IRAS 09320+6134 (see discussion below) and 
the Seyfert-2 galaxy Mrk273. An in-depth presentation of 





et al. 


(2002 


) and Ptak et al. ( 


2003) 



Compton-thick: a) the detection of an absorption turnover 
b) a flat X-ray spectrum with r ~ 1 or flatter, sug- 
gestive of a reflection-dominted spectrum, or finally c) a 
high EW (> 0.8 keV) FeKa line, indic ating a strongly ab- 



sorbed continuum (see discussion in Matt et al. 



2004 



|Akylas fe Georgantopoulosl |2009| |Gcorgant opoulos et al 



2009), for justification and explanation of these diagnos 
tics) On the basis of the above criteria, we consider 
four Compton-thick sources: IRAS 09320+6134, Mrk231, 
Arp 220, and NGC 6240. In the cases of Mrk231 and 
NGC 6240 the absorpti on turnovers have b een directly de- 
tec ted b y BeppoSAX (Braito et al. 2004 and Vignati et 



al 



1999| respectively) and therefore these can be consid 
ered as bona-fide Compton-thick sources. In the case of 
Arp 220, the spectrum is flat (r ~ 1.1) but no Fe line can 
be detected in the Chandra data. However, in the XMM- 



Newton data ( jlwasawa et al.| |2005g ) a high EW (1.9 keV) 
FeKa line at 6.7 keV has been clearly detected render- 
ing this source a highly probable Compton-thick source. 
The case of IRAS 09320+6134 (UGC 5101) is the most am- 
biguous one; we detect a flat spectrum (r ~ 1.2), al- 
beit with a large uncertainty, together with a high EW 
(~ 3.6 keV) FeKa line. The combination of a flat spec- 
trum with a high EW FeKa line rather points towards a 
Compton-thick AGN. The EW derived here is c onsistent 
with t h e EW obtained by |Ptak et al7~| ( |2003[ ). |Imanishi 
et al.| ( |2003[ ) analyse both the XMM-Newton and the 
Chandr a data finding a considerab l y lowe r EW (~ 400 eV) . 



Finally, Gonzalez-Martin et al. (2009) claim the detec 



tion of a mildly Compton-thick AGN with a column den- 
sity of N H ~ 1.4t° o4 x 10 24 cm~ 2 . Finally, in the case of 
IRAS 23365+3604, the spectrum is flat but with a large 
uncertainty. Moreover, there is no evidence for a high- 
EW FeKa line, therefore the possibility that this source 
is Compton-thick is weaker. 



2.3. Galaxies or AGN ? 

From Table 1 we see that a large number of sources (with 
the exception of Mrk273, IRAS 05189-2524 and the four 
Compton-thick AGN) are unobscured having column den- 
sities Nh < 10 22 cm~ 2 . It is possible that in a number of 
these sources the X-ray emission comes from star-forming 
processes rather than the AGN. In order to explore this 
issue, we are plotting the observed (uncorrected for obscu- 
ration) X-ray luminosity against the far-IR luminosity in 
Figure[TJ |Ranalli, Comastri fe SettI] (20031 have found a 
strong correlation between the X-ray and the far-IR lu- 
minosity for star-forming galaxies in the local Universe. 
This implies that sources with an excess of X-ray emission 
above this relation can be secure ly associated with A GN. 
The sources that fall within the Ranalli et al. (2003) re- 
lation should be either 'normal' galaxies (i.e. without an 
AGN) or highly obscured (Compton-thick) AGN. From 
Fig.[ljwe see that three sources (IRAS 05189-2524, Mrk273 
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Table 1. X-ray Properties 
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The columns are: (1) Name; A \ denotes candidate Compton-thick sources (2) Chandra Observation ID in the case where 
we performed our own spectral fits (3) Column density in units of 10 22 cm -2 (4) Power-law photon index (5) Temperature 
of the thermal component in keV (6) Energy of the FeKa line in keV (7) EW of the FeKa line in eV (8) Logarithm of the 
absorbed hard (2-10 keV) X-ray luminosity (erg s -1 ) (9) Logarithm of the fa r -infra red luminosity (40-100 fim) in units of solar 



X-ray reference: a. 


Ptak et al. 


(2003), 


Braito et al. J2004 


i f . Xia et a 


1. (2002) 


Iwasawa et al. (|2005 I x. this work. 



Kim|(|1999|) or the NASA Ex tragalactic Database ( 1 1 ) 
~?t d, |Iwe 



Ibng et al. 



Clements (2002 



and NGC 7469) present X-ray emission above that of 'nor- 
mal' galaxies. 



In order to further investigate the energy mechanism 
of the sources which lie in the window proposed by P10, 
we us e the mid-IR diagnostic diagram of |Laurent et al.| 
(2000). This traces the excess of the 5.5 //m continuum 



Table 2. AGN diagnostics 



emission with respect to the 15 /mi continuum and 6.2 /im 
PAH emission; essentially it is a tracer of hot dust from 
the AGN which has a higher 5.5 /im/15 /im flux ratio 
than the cooler dust associated with Hn regions, and low 
6.2 /xm PAH flux. The sources within the P10 box, namely 
Mrk273, IRAS 09320+6134, IRAS 10565+2448, Arp220, 
IRAS 17208-0014, IRAS 12112+0305, IRAS 23365+3604, 
show a fractional AGN contributi on of 77% 
14%, 21% and 22% respectively (I Annus et al 



, 80%, 5%, 60%, 



20071 and 



V. Charmandaris private communication) . For comparison, 
we note that the two Compton-thick AGN outside the box 
(Mrk 231 and NGC 62 40) show AGN cont ributions of 95% 
and 50% respectively ( |Armus et~alj |2007[ ). 



In Table|2] we summarize the results of the AGN diag- 
nostics used to examine the nature of the sources in the P10 
box. There is an agreement between the mid-IR and the X- 
ray diagnostics which find that 3/7 s ources host an AGN , 

pool 



if we use a 50% as a threshold for the Laurent et al 
diagnostic. We note that one so urce (IRAS 12112+ 0305) 



has been observed with the VLA (Nagar et al. 2003). This 
source, which is optically classified as a LINER, appears 
to show a radio core at 15 GHz with a radio power of 
logP « 22.6 WHz -1 . Therefore, the possibility that this 
source hosts a Compton-thick AGN cannot be securely ex- 
cluded. 



Name 


X-ray spectrum 


X-ray - IR 


mid-IR 


(1) 


(2) 


(3) 


(4) 


IRAS 09320+6134 


/ 




80% 


IRAS 10565+2448 






5% 


IRAS 12112+0305 






21% 


Mrk 273 


/ 


/ 


77% 


Arp 220 


/ 




60% 


IRAS 17208-0014 






14% 


IRAS 23365+3604 






22% 



The columns are: (1) Name (2) AGN based on a flat or 
absorbed spectrum (3) AGN based on deviation from the 



Ranalli et al. ( 2003 1 relati on (4) Contribution o f an AGN to 
the MIR according to the Laurent et al. ( 2000 1 diagnostic 



3. Discussion & Conclusions 

In Figure[2]we present the diagnostic diagram of P10. We 
see that out of the four Compton-thick sources (marked 
with crosses), two (IRAS 09320+6134 and Arp 220) fall 
within the 'high dust obscuration' window. The other 
one (NGC 6240) lies very close to the box having a 
marginally lower HCN(1-0)/GO(1-0) ratio. The fourth 
source (Mrk 231) lies away from the box showing highly ex- 
cited CO lines and a high HCN/CO J=l-0 line ratio. This 
presents no fundamental problem in the Compton-thick se- 
lection using molecular lines, as Compton-thick AGN are 
expected also to lie outside the 'window' of the diagnos- 
tic diagram of P10. Moreover, this source displays peculiar 
characteristics relative to the other three Compton-thick 
sources: although it is a bona-fide Compton-thick source, 
it is a broad absorption line (BAL) QSO. BAL is a special 
class of QSOs that show highly ionized gas flowing away 
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Fig. 1. The X-ray (2-10 keV) vs. the FIR (4-100 /xm) re- 
lation for our sources. The solid and dotted lines denote 



the best-fit relation of Ranalli et al. ( 2003 ) - and its asso- 
ciated 2(7 error, for normal galaxies in the local Universe. 
The red filled circles denote the sources which lie inside the 
box while open circles denote those outside the b ox in the 
diagnostic diagram of Papadopoulos et al. (2010). Crosses 
denote the four Compton-thick sources. 
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Fig. 2 . The diagnostic diagram of Papadopoulos et al. 
( 2010l ). High HCN(1-0)/CO(1-0) together with low CO (6- 
5)/ (3-2) ratios (region inside the box) mark the region 
where large average dust absorption likely surpresses the 
high-J (in this case J=6-5) CO emission that is emitted 
by a heavily dust-enshrouded starburst with large supplies 
of dense gas (and hence the large HCN/CO J=l-0 ra- 
tio) . Filled circles denote the excess X-ray emission sources 
(i.e. likely AGN) accordi ng Figure[T] Open ci rcles denote 
sources falling within the Ranalli et al. ( 2003 ) relation for 
normal galaxies (all four CO(6-5) / (3-2) upper limits fall 
again within the same relation). Crosses denote the four 
Compton-thick sources according to X-ray spectroscopy. 



from the central source with speeds in some cases exceed- 
ing 30000 km s" 1 . Most BAL QSOs present high obscur- 
ing columns in X-rays but still their broad-line-regions are 
viewed relatively unobscured by dust, probably meaning 
that their dense obscu ring medium is dust-free (e.g. Veron- 
Cetty & Veron 2000 ) . It is then reasonable that in the case 
of JVlrk 231 the higl7CO(6-5)/CO(3-2) ratio implies that the 
CO emission is not suppressed by dust, despite the presence 
of a large X-ray obscuring column density. 

In addition, the question that needs to be addressed 
is what is the nature of the remaining sources within 
the P10 'high dust obscuration' box. Apart from the two 
Compton-thick AGN there are five more sources presenting 
characteristically low CO(6-5)/CO(3-2) ratios in combina- 
tion with high HCN/CO J=l-0 ones. These are Mrk273, 
IRAS 17208-0014, IRAS 10565+2448, IRAS 23365+3604, 
and IRAS 12112+0305. Among them only Mrk273 has 
a known X-ray luminous, highly obscured (Njj ~ 4 x 
10 23 cm -2 ), but not Compton-thick AGN. Thus total mean 
column density of a gaseous disk around the AGN that 
would be responsible for its X-ray obscuration would be 
< N H (disk) >~ 2xN H = 8 x 10 23 cm~ 2 . On the other 
hand mean disk column densities as low as < Njj(disk) >= 
5.2 x 10 24 cm -2 are compatible with the observed suppres- 
sion of the CO 6-5/3-2 ratio in this object, which are 6.5 
times higher than the one deduced by the X-ray data. This 
discrepancy, while not very large given the totally indepen- 
dent means used to obtain < Nn(disk) > (molecular lines 
vs X-rays), remains considerable. However, we note that 
the large difficulty of measuring CO J=6-5 fluxes from the 
Earth, with a large telescope like the JCMT, will more often 
yield underestimates rather than overestimates (as point- 
ing errors become substantial with respect to narrow beams 
used) . This could tend to make molecular line-deduced and 
X-ray-deduced hydrogen column densities more compatible 
in cases like Mrk273. 

The X-ray emission of the remaining four sources in 
the box is probably dominated by star-forming processes 
according to their X-ray spectra, the Lx — L fir relation 
and the MIR diagnostics of Laurent et al. ( 2000 ) . This could 
explain why the sources do not appear as Compton-thick 
in X-ray wavelengths, as the origin of the X-rays is not 
a potentially deeply buried AGN but stellar. Thus such 
sources are not useful in deducing mean column densities 
of gaseous disks around AGN (and could be indeed found 
inside the selection "window" of P10 without invalidating 
its utility). 

As a consistency check it is interesting to discuss the Si 
9.7 /urn properties of the objects within the P10 box, as it 
is clearly impossible that these have substantial dust op- 
tical depths in short sub-mm wavelengths as suggested by 
P10, while at the same time they exhibit mo dest or low 



2006). Therefore 



Si (9.7/im) optical depths (e.g. Shi et al. 
the highly obscured objects according to P10 should present 
large Si optical depths. There are Spitzer-IRS spectroscopic 
observations for four sources in the P10 box: Mrk273, 
IRAS 09320+ 6134, Arp220 and IRAS 12112+0305 ((Annus 



et al. 



Their Si (6.7 fiva) optical depths are high: 
nffl5.47Tl)±0.3, 3.3+0.2 and 1.3+0.3 respectively. This 
clearly suggests an association of all the sources in the P10 
box with optically-thick systems at 9.7/im. 

In conclusion, our results show that the diagnostic pro- 
posed by P10 offers a potentially useful tool for the detec- 
tion of highly obscured X-ray luminous AGN . This bears 
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important implications for future facilities such as ALMA. 
However, our work identified two caveats. First, the AGN 
inside the P10 box are heavily obscured but not necessar- 
ily Compton-thick, as the case of Mrk273 implies. Second, 
there can be large contaminations by powerful star-forming 
systems. Additional X-ray spectral observations might be 
necessary in order to sift those out. More molecular line 
data, ranging from low to high-J CO lines, as well as low- 
frequency (and thus unaffected by dust exinction) density 
indicators, such as low-to-modest-J HCN and CS lines of 
nearby (U)LIRGs are needed (along with X-ray data) to 
verify (or not) this picture in local (U)LIRGs that are 
known to have X-ray luminous AGN. The now superbly 
operating Herschel Space Observatory that can d etect nu- 



merous high-J C O lines in local (U)LIRGs (e.g. van der 



Werf et al. 2010) in combination with ground based low- J 
CO and HCJN, CS observations presents an excellent op- 
portunity to do so. The large (U)LIRG sample now being 
observed as part of the Herschel Key Project HerCULES, 
in combination with the extensive ground-based low fre- 
quency CO and HCN survey will be ideal for such studies. 
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